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Abstract 



The D0 Collaboration has measured the inclusive jet cross section in pp 
collisions at y/s = 630 GeV. The results for pseudorapidities \rj\ < 0.5 are 
combined with our previous results at y/s = 1800 GeV to form a ratio of 
cross sections with smaller uncertainties than either individual measurement. 
Next-to-leading-order QCD predictions show excellent agreement with the 
measurement at 630 GeV; agreement is also satisfactory for the ratio. Specif- 
ically, despite a 10% to 15% difference in the absolute magnitude, the depen- 
dence of the ratio on jet transverse momentum is very similar for data and 
theory. 
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Typeset using REVTgX 



For reactions with large momentum transfers, quantum chromodynamics (QCD) treats 
complex proton-antiproton interactions in terms of simpler scattering processes involving 
only one constituent from each particle. Identifying these "parton" constituents with quarks 
and gluons, perturbative QCD calculates production cross sections for scattered partons 
(observed as showers or "jets" of collimated particles) that also depend on empirically- 
determined parton distribution functions (PDF) of the proton. 

This measurement compares the production rate of jets as a function of their transverse 
energy, E T , at two pp center-of-mass energies: y/s = 630 GeV and 1800 GeV. This compari- 
son reduces the systematic uncertainties and minimizes the prediction's sensitivity to choice 
of PDF. 

In the simple parton model, inclusive jet cross sections scale with y/s in the sense that 
the dimensionless quantity f(x?) = E T -E^, as a function of jet xt = does not depend 
on y/s In this model, the ratio of scaled cross sections for different energies is unity for 
all Xt- Although previous data |||3| exhibited significant deviation from this naive scaling, 
the dimensionless framework provides a useful context for comparison with QCD. The D0 
collaboration at the Fermilab Tevatron recently published the inclusive jet cross section at 
y/s = 1800 GeV using 95700 ±500 nb" 1 of data ||. This Letter presents our complementary 
measurement at y/s = 630 GeV, using a sample of 538 ± 22 nb _1 of data |J. Because the 
data at both values of y/s were collected with the same detector 0, many uncertainties in 
the results are highly correlated, and the ratio of the cross sections has greater precision 
than either of the absolute measurements. 

The differential jet cross section, d 2 a / 'dE-rdr], is measured in bins of Et and pseudora- 
pidity, 7] = — In (tan |) , where 6 is the polar angle of the jet relative to the proton beam. (In 

this formulation, the dimensionless cross section, averaged over azimuth, is ■ d 2 a / dExdrj.) 
The D0 reconstruction algorithm defines a jet by the total Et observed in calorimeter cells 
contained within a cone of radius 1Z = [(Ar/) 2 + (A</>) 2 ]2 = 0.7, where is the azimuthal 
angle. When two such clusters of cells overlap, they are merged into a single jet if they share 
more than 50% of the E? of the lower- Et cluster; otherwise, they are split into two separate 
jets, each defined by its own rj-cj) centroid and Et value @]. 

The online trigger requires at least one jet above a set threshold. The offline data selection 
procedure, which suppresses backgrounds from electrons, photons, noise, and cosmic rays, 
follows the methods used in the 1800 GeV analysis [@||. The efficiency of jet selection is 
approximately 96% and is nearly independent of jet Et- To maintain precision in jet Et-, 
a vertex requirement removes jets resulting from pp interactions more than 50 cm from the 
center of the detector, thereby reducing the total efficiency to 82%. The uncertainty on the 
cross section associated with all efficiencies is < 0.5% ||. 

Jet energies are corrected jL(J for the energy response of the D0 calorimeter to hadrons, 
the broadening of the hadronic shower, and energy from multiple interactions, calorimeter 
noise, and the underlying event (fragmentation of the spectator partons). The response 
correction increases the Et of jets by 22% for measured calorimeter Et of 20 GeV, and by 
15% for jet Et above 100 GeV. The 1% showering correction recovers the net energy lost 
when hadrons from inside the R=0.7 cone deposit energy outside it as they interact within 
the calorimeter. Calorimeter noise, from electronics and from uranium activity, contributes 
on average 1.6 GeV of Et to each jet. The underlying event contributes 0.6 GeV of Et 
to each jet at y/s = 630 GeV, compared to 0.9 GeV at y/s = 1800 GeV. The corrections 
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offset one another, so that a jet's measured Et typically increases by 12% to 14% after 
implementing all energy scale corrections. Uncertainties in the corrections for noise and 
response dominate the systematic uncertainty of the final result. 

Both detector imperfections and random fluctuations in shower development of individual 
jets within the calorimeter result in the smearing of a jet's Et about its true value. The 
finite Et resolution shifts the observed cross section to higher Et, especially in the most 
steeply falling regions of the distribution. The measurement of jet resolution as a function 
of Et and the unsmearing procedure follow the steps described in Ref. Q. The unsmearing 
correction is larger at 630 GeV than at 1800 GeV because the cross section is significantly 
steeper at the E T values of interest. 

Figure |l| depicts the inclusive jet cross section at y/s = 630 GeV in the pseudorapidity bin 
|r/| < 0.5. Each data point indicates the Et at which the cross section within that bin has its 
average value. The bin widths are chosen to match the bins in xt from the y/s = 1800 GeV 
analysis. Table | reports the bin ranges, point positions, and uncertainties. The solid 
line in Fig. [I] indicates the result of a calculation using the jetrad next-to-leading-order 
(NLO) partonic event generator ||TT| and the CTEQ3M PDFs [|12|]. The renormalization and 
factorization scales are set to /i = i?™ ax /2, where i?™ ax corresponds to the Et of the leading 
jet in an event. 
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FIG. 1. The inclusive jet cross section at y/s = 630 GeV, integrated over azimuth and averaged 
over 1 77 1 < 0.5. The shaded band corresponds to the systematic uncertainties in the measured cross 
section and the solid line shows a prediction from NLO QCD. 



Figure § compares the cross section to the NLO QCD prediction in greater detail. The 
"baseline" renormalization and factorization scales are set to fi = E T ia,x /2; additional lines 
in Fig. |2| indicate the predictions that result from changes in either PDF or \i relative to 
the baseline prediction specified for that pane. The shaded regions in Fig. |2| indicate the 
one standard deviation systematic uncertainty of the measurement, and the vertical bars 
indicate the statistical uncertainty. The first prediction, generated with the MRST BT3 
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PDF, is shown to best reproduce the absolute magnitude of the data, but the CTEQ4HJ 
14] curve in the second pane appears to provide the closest match in shape. Changing 
(j, modifies both the normalization and the shape of the predictions, as seen in the third 
pane. We quantify the agreement between the data and the various predictions with a x 2 
comparison, as described below. 
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FIG. 2. The inclusive jet cross section at y/s = 630 GeV compared to several NLO QCD 
predictions. Error bars indicate statistical uncertainties and shaded bands correspond to systematic 
uncertainties. The horizontal lines at zero indicate the baseline prediction that is named in each 
pane; additional lines indicate theoretical variations relative to the baseline. 

Combining the results from this Letter with those of Ref. Fig. ||] displays the ratio 
of dimensionless jet cross sections as a function of xt- The observed ratio ranges from 1.48 
to 1.85, depending on the value of xt- The largest uncertainties arise from the corrections 
for response and noise, and the rest primarily from resolution and luminosity. Although 
the systematic errors on the individual measurements range from 10% to as much as 30%, 
strong correlations reduce the uncertainty on the ratio to values as small as ±5.4%. The 
two final columns of Table | provide the numerical results for the ratio. 

As shown in Fig. |3], NLO QCD predictions for the ratio lie systematically above the 
data throughout most of the measured xt range, in particular between xt of 0.1 and 0.2, 
where the ratio has the smallest statistical uncertainty. Choice of PDF has little effect 
on the prediction — only the renormalization/factorization scales change the prediction 
appreciably. 

A covariance matrix x 2 comparing data and theory provides a measure of the probability 
that the theory describes the observed results. To verify that our covariance matrix, built 
mostly from correlated systematic uncertainties, produces results that are consistent with 
a standard x 2 distribution with 20 degrees of freedom, we generated an ensemble of 20 
million experiments using a Monte Carlo program. Each statistical and systematic error was 
simulated and varied randomly using appropriate correlations in and y/s. Systematic 
errors were not necessarily assumed to be Gaussian distributed; some numbers were drawn 
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TABLE I. Inclusive jet cross section at yfs = 630 GeV and the ratio of dimensionless cross 
sections / 630 (x T )// 1800 (x T ), where / (x T ) = E T ■ and x T = 2E T /^. The cross sections are 

all integrated over azimuth and averaged in the range \rj\ < 0.5. 
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FIG. 3. Ratio of dimensionless jet cross sections (numerator ^/s = 630 GeV, denominator 
t/s = 1800 GeV) compared to NLO QCD as given by jetrad. Error bars indicate statistical 
uncertainties and shaded bands correspond to systematic uncertainties. 



from uniform probability distributions, as appropriate. The x 2 comparisons (between the 
original input and each of the final, randomly- varied distributions) is in excellent agreement 
with the shape of the x 2 function for 20 degrees of freedom. We find that the standard 
probability obtained from an integral of the x 2 distribution provides an appropriate vehicle 
for comparing data with predictions. 

Table H reports both the x 2 values and the x 2 probabilities for the comparison of the 
data with different NLO QCD predictions. The inclusive jet cross section at \/s = 630 GeV 
is consistent with all the tested PDFs and fi scales, with but two exceptions. For the ratio 
of cross sections, there is no significant difference in shape between data and theory, and 
essentially all predictions lie within an acceptable range. The overall results in Table |TJ 
indicate reasonable agreement between the ratio and NLO QCD. 

We performed a second test to quantify the observed difference in the absolute magni- 
tudes of the predicted and observed ratios, without particular regard to the shapes of the 
distributions. Using the covariance matrix and assuming that the value of the ratio is a 
constant with respect to xt, we found the best-fit horizontal line for the data. The x 2 value 
that results from a comparison of this single point to the equivalently-calculated theory 
point yields the probabilities listed in the final column of Table [IT]. In every case, discarding 
the information on shape in favor of a comparison of absolute magnitude results in poorer 
agreement between data and theory, particularly for the often-favored scale of fi = i?™ ax /2. 

In conclusion, we have measured the inclusive jet cross section at two center-of-mass 
energies, 630 GeV and 1800 GeV. Both the published data at 1800 GeV gj] and the data 
presented here at 630 GeV are generally well-described by NLO QCD, with the exception of 
predictions using CTEQ3M(2i?™ ax ) and MRSTGD PDFs. In the ratio of dimensionless cross 
sections at the two energies, experimental uncertainties are much smaller and differences 
in the predictions from choice of PDF are less important. NLO predictions for the ratio 
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PDF 


A* 


630 GeV C. Sec. 


Ratio 


Norm. 






x 2 


Prob. 


x 2 


Prob. 


2 

A. 


Prob. 




2 . ^max 


40.5 


0.43% 


17.9 


59.4% 


3.33 


6.81% 




Ex 


25.9 


16.8% 


21.6 


36.2% 


7.13 


0.76% 


CTEQ3M 


E™ x /2 


30.4 


6.37% 


20.5 


42.5% 


9.56 


0.20% 




E™ x /4 


27.5 


12.2% 


15.1 


77.4% 


1.45 


22.93% 


CTEQ4M 




24.1 


23.8% 


22.4 


31.9% 


10.67 


0.11% 


CTEQ4HJ 


E™ x /2 


18.9 


52.5% 


21.0 


40.0% 


13.21 


0.03% 


MRST 


£™ ax /2 


22.6 


30.7% 


22.2 


33.0% 


12.60 


0.04% 


MRSTGU 


Ef ax /2 


14.9 


78.2% 


19.5 


48.7% 


11.07 


0.09% 


MRSTGD 


Ef ax /2 


51.8 


0.012% 


24.1 


23.9% 


12.92 


0.03% 



TABLE II. x 2 comparisons for the cross section at y/s = 630 GeV (20 degrees of freedom), the 
ratio of cross sections (20 degrees of freedom), and a comparison for the ratio involving only the 
absolute magnitude (one degree of freedom). 



exhibit satisfactory agreement with the shape of the observed ratio. In terms of only the 
magnitude however, the absolute values of the predictions lie significantly higher than the 
data, especially for the standard scale fi = i?™ ax /2. 
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